General Theory of Lee- Yang Zeros in Models with First- Order Phase Transitions 
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We present a general, rigorous theory of Lee- Yang zeros for models with first-order phase transitions that 
admit convergent contour expansions. We derive formulas for the positions and the density of the zeros. In 
particular, we show that for models without symmetry, the curves on which the zeros lie are generically not 
circles, and can have topologically nontrivial features, such as bifurcation. Our results are illustrated in three 
models in a complex field: the low-temperature Ising and Blume-Capel models, and the g-state Potts model for 
q large enough. 
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Almost half a century ago, in two classic papers |]J], Lee 
and Yang studied the zeros of the Ising partition function in 
the complex magnetic field plane, showed rigorously that the 
zeros lie on the unit circle, and proposed a program to ana- 
lyze phase transitions in terms of these zeros. A decade later, 
Fisher |2] extended the study of the Ising partition function 
zeros to the complex temperature plane. Since that time, there 
have been numerous studies, both exact and numerical, of 
the Lee-Yang and Fisher zeros in a wide variety of models 
01 . However, with a few notable exceptions [4], remarkably 
little progress has been made in extending the rigorous Lee- 
Yang program. This is due to the fact that rigorous statistical 
mechanics has relied almost exclusively on probabilistic tech- 
niques which fail in a complex parameter space. In this Letter, 
we adapt complex extensions I.SatZD of Pirogov-Sinai theory 
|0] to realize the Lee- Yang program in a general class of mod- 
els with first-order phase transitions. 

The purpose of this work is threefold. First, it is of in- 
terest to establish the mathematical foundation of a program 
that has been so central to statistical physics. Second, our 
theory gives a novel physical interpretation of the existence 
and position of partition function zeros by relating them to the 
phase coexistence lines in the complex plane. Finally, from a 
practical viewpoint, our theory provides a framework for the 
interpretation of numerical data by allowing explicit, rigor- 
ous computation of the position of the zeros. Indeed, we find 
rigorous results which clarify many ambiguities in published 
data. Specifically, in models without an underlying symme- 
try, we prove that the zeros generically do not lie on circles, 
even in the thermodynamic limit. This applies, in particular, 
to the Blume-Capel and Potts models in complex magnetic 
fields; see |0 1 1 Oil for heuristic studies of these models. We 
also prove that the curves defined by the asymptotic positions 
of the zeros can have topologically nontrivial features, such 
as bifurcation and coalescence, and show that these features 
correspond to triple (or higher) points in the complex phase 
diagram. 

The results to be stated next are rather technical. Roughly 
speaking, they say that, for models with a convergent contour 
expansion, the partition function can be written in the form 
([0 and its zeros are given as solutions to equations Q and 



(|3}. Readers unfamiliar with rigorous expansion techniques 
are encouraged to see the concrete examples following the 
main results. 

Main Result: Consider a d-dimensional lattice model with 
n equilibrium phases whose interaction depends on a com- 
plex parameter z. Suppose d > 2 and that z is in the region 
(typically, a large disc or the entire C) where the model ad- 
mits a contour representation with strongly suppressed con- 
tour weights. Under suitable conditions |5[ |fj |8|], there are 
complex functions fi = ft(z), £ — 1, . . . , n, such that the 



partition function in a periodic volume V 
temperature (3 can be written as 



L at inverse 



-0ftV 
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Here Lrj is of the order of the correlation length, / = 
minfc 'Stcfk, and qg is the degeneracy of the phase £. Physi- 
cally, fi can be interpreted as metastable free energies with 
the stability of the £ th phase being characterized by the condi- 
tion diefg = f. If £ is stable, fg is just the free energy of the 
system with boundary condition £ 111 ill . In the region where 
£ is not stable, fg is constructed as a smooth extension of fg 
from the stable region. Clearly, fg depends on the parameters 
of the model, but not on L. 

Eq. Q can be used to locate the zeros of Z£ er analytically. 
Excluding a neighborhood of size 5 L ~L-( d -V of the triple 
or higher coexistence points lll2ll and assuming a degener- 
acy removing condition fl3ll . each zero of Z^ T lies within 
0(e~ L / L °) of a solution to the equations 



< Re/ft for all k ± £, m, (2) 



(3V(^smfg — 3m/ m ) = n mod 2-7T 



(3) 



for some £ ^ m, where f£ s L — ft — ((3V) 1 logq^. In fact, 
the solutions to (|2ji and Q and the zeros of Z£ er are in one- 
to-one correspondence. As a consequence, the zeros of Z^ r 
asymptotically concentrate on the phase coexistence curves 
Uzfg = Re/ m with the density ±(3V\{d/dz){fg - f m )\. In- 
side the (^-neighborhood of the multiple coexistence points, 
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both the analysis and the resulting equations for the zeros are 
more complicated; see |7|] for details. However, it turns out 
that all but a uniformly bounded number of zeros (out of the 
total of order L d ) can be accounted for by the simple equa- 
tions (0 and (0. 

The proof, which appears elsewhere 0], is technically com- 
plicated, but the main idea is simple. The key input is a com- 
plex version of methods developed mainly in the context of 
finite-size scaling lEH^Bl leading to equation ([0 and similar 
expressions for the derivatives of Z^ el . Equations and (0 
for the zeros of Z£ el arise from "destructive interference" of 
two terms, qie~P^ v and q m e~^^ mV , in the sum in Q. Out- 
side the 8l -neighborhood of multiple coexistence points, all 
other terms are negligible. 

To illustrate our result, we will discuss three specific mod- 
els in the presence of a complex external field. 

Ising Model: The nearest-neighbor Hamiltonian is 

0H = - J } <7 x <jy - h} a x . 

Here a x G {— 1,+1}, the coupling J > is taken 
large enough to ensure absolute convergence of the low- 
temperature expansion, and h is the complex external field. 
Neglecting the error term, Eq. Q becomes 



This leads to the following equations for the zeros: 
»e(/+-/_) = 



9m(/+-/-) 



(2fc — 1)71 

0V 



,.,V. 



(4) 



(5) 



Inserting the low-temperature expansions of f±, 

0f± = Th-dJ- e - 4dJ e T2h + R{±h), 



where R{h) and its derivative are both 0(e 

e ie ", with 



-4(2d-l)J' 



, we 



find that the zeros occur at e 



2h 



9 k = 



(2fc-l)7 

V 



2e 



-4dJ 



3in ((^_ll£) + 0( 



(k e -4(2d-l)J 



k = 1, . . . , V. Moreover, the h <-> — h symmetry can be 
used to prove that condition (0J is equivalent to 5Re h = 0, 
guaranteeing that the zeros of Z^ 1 lie within an 0{e~ L / L °)- 
neighborhood of the unit circle. The h <-> —h symmetry of 
the partition function then allows us to conclude that, for large 
L, the zeros lie exactly on the unit circle; see the end of the 
Blume-Capel section for details of an analogous argument. 
This gives an alternative proof of the Lee-Yang circle theo- 
rem at low temperatures |7]. We stress that symmetry is the 
key factor here; in the absence of symmetry, (|4} does not in 
general lead to circles. 

Blume-Capel Model: The Hamiltonian fl4ll is 



(3H 



J(°x - Oyf - J2(^crl + ha x ) 
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FIG. 1: The 128 zeros of the partition function of the Blume-Capel 
model in the complex e h plane, for the 8x8 periodic square grid 
at e _4J = 1/16 and e A = (a) 0.9, (b) 0.94, (c) 1, (d) 1.07. The 
actual zeros lie within a distance of order e~ lt)J +e~ L / L ° of those 
depicted. For A < A^T ~ — e~ 4J , the outer region of the + phase 
is separated from the inner — phase by an annular region of the 
phase (a); asymptotically, the zeros lie on the boundaries of these 
regions. As A increases through A~, the two boundaries coalesce on 
the left hand side, leading to bifurcation for A > A^T . The common 
boundary grows (c,d) and, eventually, at A = A^ &s e _4J , the 
phase disappears and bifurcation terminates. For A > Aj, all zeros 
lie on the unit circle. 



with spins a 6 { — 1, 0, +1}, real parameters J > 0andA,and 
complex field h. For J large, the real (A, /i)-phase diagram 
features three phases labeled by +, 0, and — , each with an 
abundance of the corresponding spin. 

The zeros of this model are shown in Fig. 1. Note that the 
zeros have a non-uniform distribution, forming curves of non- 
circular shape, and that for A in a certain interval (X~ , A 4 "), bi- 
furcation (i.e., splitting of the curve) occurs. In the remainder 
of this section, we rigorously establish these features for large 
J. Before beginning our analysis, we remark that in 1 9] a phe- 
nomenological theory of partition function zeros based on |6] 
was developed and then applied to the Blume-Capel model. 
In contrast to our approach, that of |9] gives no quantitative 
estimate of approximations or errors, and it misses certain im- 
portant qualitative features, namely the bifurcation. 

Fix J large and let e h = z = re 10 . Our analysis is done 
in two steps. First, we focus on the unit circle, r = 1, and 
identify Xf and the position e ±ie d x ) of the splitting points. 
Then we extend the analysis to all r. 

The J = oo phase diagram has three ground states, a 
1 , 0, 1, with energy densities —A + h, 0, —A — h. The large- 
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J expansions of the free energies are 

e -^/ + = ze A exp {I e -A-4,7 + 2 ^ e ~2X-6J + 
fi -/S/_ = i e A exp { ze -A-4,7 + 2z 2 e -2A-6J + j 

e-^=exp{(z+i)e A - 4 ' / + 2(z 2 + ^) e 2 ^ w + i? (z)}, 

which, for brevity, we write only for d = 2 11511 . Here i?± (z) 
and Ro(z) and their first two derivatives are all 0(e~ 8J ). 
There are no phase degeneracies, qo = q± = 1. 

By Eq. 0, the analysis of the loci of the zeros requires 
comparison of 5Re/+, 5ft e/_, and Sfte/o. On the unit circle, 
5Re/+ = 5fte/_. Thus it suffices to study the sign of A± = 
SRe/?/± - 3?e/3/ . For |A| > e" 4 ' 7 , sgn(A±) = -sgnA. So 
let A = 0(e- iJ ). Then, for J large, 

A± = -A + e- iJ {2e x -e- x ) cos(9 + ©(e^ 67 ), (6) 

so that (d/d0)A± < for 9 G (0,tt) JH and, similarly, 
(d/dX)A± = — 1 + 0(e~ 4J ). This implies the existence of 
A± = ±e" 4J + C(e- 6J ) and C (A), with C (A) G (0,tt) for 
all A G (Aj, A+), such that, on the unit circle, is the only 
stable phase for all 9 when A < A~ and for \9\ < 9 C (X) when 
A G [A~,A+], whereas ± are the only stable phases in the 
complementary region of (A, 9). Moreover, 9 C (X) decreases 
with A and 9 C (X) — ► (resp., it) when A | A+ (resp., A j Aj). 

Now let r be arbitrary. We have (d/dr)§tef± (z) = ±r~ 1 + 
C(e" 4J ) and (d/dr)5Re/ (z) = ©(e" 4 - 7 ). Using the symme- 
tries of the model, 

SRe/±(z) =5Re/ T (z- 1 ) and Re/ (z) = aSe/o^ 1 ), 

it follows that there is a function i— ► r(0), < 1 — r{9) < 
0(e- iJ ), r{9) = r(-9), such that - is stable for r < r(0), 

is stable for r(9) < r < l/r(9), and + is stable for r > 
l/r{9). Notice that r{9) = 1 for \9\ > 9 C {X) when X~ < X < 
A+ and for all 9 when A > A+. 

Consider now L ^> Lq and suppose there is a partition 
function zero at zq — pe 1 ^ . If r(i/j) < 1, then the zero lies 
close to one of the curves defined by the equations A± = 0. 
We claim that these curves are non-circular and that the zeros 
do not maintain a uniform spacing along them. Indeed, set 
A = for simplicity and observe that 

e A± = |zT 1 exp{z ±1 e- 4J + 2z ±2 e- 6J + 0(e- 8J )}|. 

Replacing z^ 1 by x + iy, the equation A± = and the ex- 
pansion of the exponential up to (D(e~ 8J ) yield 

x 2 + y 2 = 1 + 2xe~ 4J + 4(x 2 - y 2 )e- 6J + <D(e~ 8J ). 

This is an ellipse centered at e~ 4J + 0(e~ 8J ) with semiaxes 

1 ± 2e~ 6J + 0(e~ SJ ). To determine the density of zeros, we 
compute \(d/dz)(f± — fo)\ and easily verify that it is non- 
constant on the above ellipse. 

If, on the other hand, r(tp) = 1 IHIL then, for L large 
enough, the zero necessarily lies exactly on the unit circle, 
p = 1. Indeed, by (0, 0, and the degeneracy removing con- 
dition 1 13;], the distance between two adjacent zeros is of order 




0~ 1 -1 ~~ (T~ 



FIG. 2: Complex-e' 1 dia gram showing the zeros of the Potts model in 
a three-dimensional periodic box of size V = 1000 with parameters 
q = 25, e 3J /q = (a) 1.185 and (b) 1.155. In each case, there are 
1000 zeros distributed on three non-circular arcs, labeled A, B, and 
C, with those on A and B denser than those on C. The outer region 
corresponds to the ordered magnetized phase, while the regions left, 
resp., right, of arc C contain the ordered non-magnetized and disor- 
dered phase. For V large, arc C shows up first at J — J', it passes 
through zero at J = J^" 1 , and it disappears at J ~ Jc~ 2 . 

L~ d . But we also have \p - r(ip)\ = \p - 1| < O(e- L / Lo ), 
and if p ^ 1, then by symmetries of the model, there would 
be another zero at z^ 1 = p~ x e 1 ^ . However, |z — z$ 1 \ < 
O(e _i/io ) <C L~ d , a contradiction. A similar argument 
proves a "local" version of the Lee-Yang theorem 11811 in a 
large class of models for which the standard, "global" theo- 
rem fails. 

Potts Model: The Hamiltonian is 

/3H = -J^2S !Tx!!Ty -hJ2^,i- ( 7 ) 

with spins a G {1, 2, . . . , q}, real coupling J > 0, and com- 
plex field h. For h = this is the standard Potts model, with 
a q-fold degenerate ordered phase at large J and a disordered 
phase at small J, coexisting at « ^ log q. The transition is 
first-order for large q QjJ, while it is presumably second order 
for q < q c (d). For h ^ and q large, the phase diagram was 
determined first by formal expansion 1 20 ] , and recently by rig- 
orous probabilistic methods |21]. The Lee- Yang zeros of @ 
were studied numerically in 1 10], where it was suggested that 
the zeros lie on almost circular curves slightly outside the unit 
circle, for J both above and below J*. While, by three-phase 
coexistence, this turns out to be incorrect for J < Jjf (see 
Fig. 2), we prove that this is indeed the case for J > J%, thus 
resolving a controversy in fioll . 

The model Q has three phases: the disordered phase (D) 
with degeneracy qjj = 1, and two ordered phases: a magne- 
tized (M) and a non-magnetized (O) phase, with degenera- 
cies gjvr = 1 and qo = q — 1, characterized by abundances of 
cr z = 1 and a z = const. ^ 1, respectively. Let us abbreviate 
z = e'\ n d = d(2d - 1), Q ( z k) = q - 1 + z\ and Q z = Q ( z 1] . 
The free energies are given 12111 by 
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e-^ = Q z cxp{d(e J - 1)Q ( V/Q 2 Z + n d {e J - l) 2 Q z ^/Q 3 z - {n d + l/2)(e J - l) 2 [Q^/Q 2 z } 2 + <D(l/q 3 - 4 / d )} 
e -f3f M = zexp { dJ + e -^J {QzZ -i + de -( id -V J (Q 2 z + e J Q^)z- 2 - (d + l/2)e~ AdJ Q 2 z z~ 2 + 0{l/q 3 - 2 ' d )} 
e-^o = exp{dJ + e- 2dJ (Q z - 1) + de"^" 1 ' J (Q 2 Z + e J Q^) - (d + \/2) e - idJ Q 2 Z + G(l/ g 3 - 2 / d )}. 



The zeros of the periodic Potts partition function are de- 
picted in Fig. 2. In particular, for J_ < J < J + (where 
J + = J% and J_ « Jq~ 2 ), the loci do not lie on a sin- 
gle closed curve but rather split the complex plane into three 
pieces, corresponding to the regions of stability of the three 
phases above. The number of zeros on the inner arc is roughly 
V I (2Trq), so one needs to take V quite large and tune J to fall 
inside the narrow window ( J_ , J + ) to find any interior zeros. 
This explains why these zeros were not detected in previous 
numerical work 1 10]. 

Despite their appearance, none of the curves in Fig. 2 is 
a circle. This is verified by finding the coexistence curves (0 
for three distinct pairs k, £ G {D, M, O}. When J > J«, only 
the phases M and O are relevant, and the asymptotic location 
of the zeros is given by Sfte/o = 5ie/A/- For z = re 10 , this 
easily implies 

r = 1 + qe- 2dJ (l - cosfl) + 0(l/q 2 ), 

so that for J > J« and V — > oo, all zeros with |0| > 
are asymptotically outside the unit circle. By invoking argu- 
ments similar to 1 16], this extends to all 9 12111 . There are two 
finite-volume corrections: an outward shift of order 1/V due 
to fo > Pol ( see Eq. and an error 0(e~ L l La ) coming 
from Q. Since 1/V > 0(e~ L/L °), this proves the initial 
numerical observation in II 1011 . 

To make the interesting features clearly visible, Figs. 1 
and 2 were drawn for values of e~ 4J and q for which we have 
not proved convergence of our expansions. However, as es- 
tablished above, all the depicted behaviors indeed occur once 
e~ J (or 1/q) and e ~ L / L ° are small enough. 

In summary, we identify the loci of complex zeros with the 
complex phase coexistence curves. For particular models, we 
use this identification to map the precise location of these ze- 
ros. We find that, in general, the loci are non-uniform and that 
the resulting curves are non-circular; if more that two phases 
are present, the curves also have bifurcation (i.e., splitting) 
points. 

The research of R.K. was partly supported by the grant 
GACR 201/00/1 149. 
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